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Abstract 

An effective Lagrangian for the light quark in the field of a static source is derived 
systematically using the exact field correlator expansion. The lowest Gaussian term 
is bosonized using nonlocal colorless bosonic fields and a general structure of effective 
chiral Lagrangian is obtained containing all set of fields. The new and crucial result 
■ is that the condensation of scalar isoscalar field which is a usual onset of chiral 



symmetry breaking and is constant in space-time, assumes here the form of the 



confining string and contributes to the confining potential, while the rest bosonic 



fields describe mesons with the qq quark structure and pseudoscalars play the role of 



D . Nambu-Goldstone fields. Using derivative expansion the effective chiral Lagrangian 

' is deduced containing both confinement and chiral effects for heavy-light mesons. 

The pseudovector quark coupling constant is computed to be exactly unity in the 
^ . local limit, in agreement with earlier large Nc arguments. 

1 Introduction 

It was understood long ago that chiral symmetry breaking in QCD is responsible for the 
low mass of pions and therefore the low-energy limit of QCD can be adequately described 
by the effective chiral Lagrangians 0. 

In this approach the Nambu-Goldstone particles are described by local field variables 
and the resulting Effective Chiral Lagrangian (ECL) is local. 

The most general and practically useful form of chiral Lagrangian was given in and 
contains around ten phenomenological parameters (14 to the 4th order in p), to be found 
from experiment. 

Being successful in describing low-energy processes with Nabmu-Goldstone mesons, 
ECL has two major defects: Firstly, it does not take into account the quark structure of 
mesons, and consequently e.g. the formfactor computed in ECL can display in the meson 
only mesonic degrees of freedom. 
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Secondly, ECL completely disregards the phenomenon of confinement, which is also 
important at small energies, and hence degrees of freedom of vacuum gluons, creating 
confining string are not taken into account. 

There have been attempts to cure the first defect, namely the model Lagrangians have 
been suggested which take into account both quarks and mesons 0, 3- In particular the 
instanton model of the QCD vacuum has been used to derive the ECL and the Quark- 
Meson Lagrangian(QML) [§. As a result interesting interconnections of quark and chiral 
degrees of freedom have been demonstrated in the example of the nucleon 

However instantons are suppressed in the reahstic vacuum of QCD H, || and more- 
over the internal consistency of the instanton vacuum without confinement is seriously 
questioned |TT[ . 



Moreover instantonic vacuum lacks confinement and therefore cannot cure the second 
defect of ECL and QML, therefore model ECL and QML obtained in |^ disregard 
confinement fully. 

Finally, the Effective Quark Lagrangian (EQL) obtained after averaging over gluon 
degrees of freedom, contains in principle infinite number of terms with growing number of 
quark fields. It was shown in that specifically for instanton vacuum all higher terms 
are of the same order, while only the lowest term (so-called 'tHooft Lagrangian) is taken 
in the standard instanton lore. 

It is the purpose of the present paper to start a new and systematic approach to the 
derivation of QML and ECL from the first principles - the QCD Lagrangian. In doing so 
the most important is to keep gauge invariance at every stage, therefore we shall consider 
the simplest gauge-invariant system of a light guark in the field of a static antiquark 
(generalization to only light quark systems will be done later). 

After averaging over gluon fields one obtains EQL with infinite number of terms, 
containing as kernels irreducible gluon Field Correlators (FC). Recent measurements for 
the realistic QCD vacuum have shown that FC create a hierarchy, where the lowest 
(Gaussian) correlator is dominating while the next (4th order) correlator contributes 
around 1% to the static QQ interaction (note that this situation is drastically different 
from that of instantonic vacuum, where higher correlators are equally important). We 
assume that a similar hierarchy should be present also in our problem of a light-heavy 
quark, which allows us to study the resulting EQL term by term, paying most attention 
to the lowest, 4g-piece. The next step is the bosonization procedure, i.e. an identical 
transformation introducing nonlocal colorless bosonic fields, having different Lorentz and 
flavour indices due to the use of Fierz transformation. 

A special attention deserves the 6q term, where bosonization may be done introducing 
bosonic fields for qq combination or else introducing baryonic fields for 3q combination. 

As a result one obtains QML or Quark-Baryon Lagrangian (QBL) in the most general 
and rigorous form. These questions will be considered in a separate publication |jl3 . 



As the next step one can use the stationary point analysis to obtain nonlinear equations 
for effective bosonic fields. 

This is done in a rigorous way and the resulting equations contain for the scalar- 
isoscalar part the same equation as was derived previously fl^ in a gauge-invariant Dyson- 
Schwinger approach not using bosonization. The careful study of that equation in WM- 



16| has shown that it describes in the light-heavy system both confinement and chriral 



symmetry breaking, which coexist. This fact means that confinement in the language of 
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effective meson fields, enters in the form of condensate of the scalar field inside the string, 
while other fields describe subdominant features of the qQ dynamics. 

The paper is organized as follows. In the next section the averaging over vacuum 
gluonic fields is done and EQL is obtained. Special attention is devoted to the gauge 
invariance and parallel transporters necessary to ensure it for nonlocal qq combinations. 
In section 3 the bosonization procedure is done and Fierz transformation is introduced to 
obtain final QML with proper classification in Lorentz and flavour indices. 

In section 4 the resulting ECL is obtained and the stationary point equations are 
derived and compared to the previously obtained in the Dyson-Schwinger approach. In 
section 5 the derivative expansion of the nonlocal bosonized Lagrangian is done and non- 
local forms of the lowest 2nd and 4th order terms in this expansion are obtained. Keeping 
only pion field in addition to quarks, one derives in section 6 the effective Lagrangian 
which appears to have the expected form with the pseudovector quark-pion coupling con- 
stant g\. The latter is equal exactly to one, in agreement with earlier large argument, 
when local approximation is made. The concluding section is devoted to the discussion 
of confinement and chiral properties of the resulting ECL in the heavy-light meson case. 

2 The Effective Quark Lagrangian 

Consider the QCD partition function in the Euclidean space-time 

Z = j DAD^D^+e-^'^^^^+I ^^-^(id+i^i+gA) f-^d^x (^-^^ 

where So{A) = ^ m is the current quark mass, and the quark operator 
■^i^aaix) has flavor index /(/ = 1, ...nf), color index a(a = 1, ...Nc) and Lorenz bispinor 
index a{a = 1, 2, 3, 4). 

The next step is to integrate over DA^^ with the weight So{A), this is the gluon vacuum 
averaging which is denoted by ( )a- Before doing this however one should choose the 
gauge-invariant system and using an appropriate gauge, express A^ through the field- 
strength operator F^^, which would finally appear in gauge-invariant combinations - Field 
Correlators (PC) (see for review and more discussion), namely 

))a ^ A(") (2) 

where parallel transporters are defined as 

^c{x,y) = Pexp{ig A^,dz^) (3) 

JC{x,y) 

and the open contour C (x, y) connects points x and y and can be arbitrary otherwise. 

To achieve this goal one can use the so-called contour gauge ^ which is especially 
convenient in the case of light quark moving in the field of a static antiquark. One has 
for the contour z^{s^ x) starting at point x and ending dXY = 2;(0, x) 

A,{x) = £ ds^^^^^^j^F^Ms)) ^ £ dT,^,{z)F^,{z). (4) 
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In particular case, when the contour z^{s,x) goes along the shortest (straight) way to 
the X4- axis and then along this axis to some point Y, which can be at - 00, one has the 
so-called modified Fock-Schwinger gauge ||20|, which was extensivelly used in 0] to get 
EQL. Here one has 



a^{u)duiFi^{u,Xi) 



(5) 



and a^iu) = 1 while for /x = 1, 2, 3, is equal to 



aiu) 



Ui 



Xi 



1,2,3. 



It is convenient to write all expressions in a gauge-invariant way, using the property ||19||, 
that given by (H) is identically equal to unity when the contour C{x, y) lies on z^{s, x) 
or z^{s,y). Therefore one can define gauge-covariant operators referred to the point F, 

ip^^\x) = <l>cix,Y)^{x) = ^Y,x)^{x) 

=^+(x)<l>(x,F) 

F^l\x) = <^(Y,x)F^,{xMx,Y). (6) 

Here the contour C{x,Y) in (j){Y,x) goes along z{s,X) from Y to x, and in opposite 
direction in y). 

For the field correlators referred to the same point Y one can write {ab, cd- fundamental 
color indices, Lorentz indices are suppressed for simplicity reasons) 



g\{F^''\x)U{F''^\y)U 



^-^9'{tr{F^''\x)F^'^\y))) 



(7) 



^7^((F(^)(x)),,(F^(y)),,(F(^)(^))e/) = 6,J,Ja) 



-g^tr{xyz) 



+ SbeSfc^ad 



g^tr{xyz) 



+ 



gHr{xzy) 

>c(iVc'-l)' 

NcgHr{xzy) 

m-1 



where notation is introduced e.g. tr{xyz) = {tr{F^'^\x)F'^^\y)F^'^\z))) . 
Derivation of ®,(|D is given in Appendix 1. 

Let us first concentrate on the bilocal correlator (|^). From (^, (j^) it is clear that in the 
average value of {A^{x)Ay{ii)) the arguments of F{z{s, x)) and F{z{s, y)) are separated by 
the distance r ~ T^, where Tg is the gluonic correlation length [^], which was measured 
on the lattice |2^ and estimated analytically |]23[ to be Tg ~ 0.2fm (or even smaller, 
if data on gluelump masses |^ are used). For such distances r ~ T^, which satisfy 
r <^\x — Y\, \y — Y\ or for the gauge (13), r ^ |x|, |y|, one has an estimate 



(trF(^)(x)F(^)(y)) = {tr{F{x)^{x,y)F{y)^{y,x))) + 0{-^J- 



X^ y2 ' 



(9) 



where the correlator on the r.h.s. of (^ is connected by straight lines from x io y. A 
similar estimate holds for the triple FC (^, and in what follows we shall use the straight- 
line form (P|) which is independent on the position of the reference point Y . 
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For that correlator one can use the general representation found in P3 

^^{trF^,{x)^{x,y)FpM^{y,x)) = D{x - y){5^p5,„ - + ^Sp. (10) 

where A*^^) has the structure of a full derivative and therefore does not contribute to 
confinement, its nonperturbative (NP) part is much smaller than that of D and we shall 
omit it. 

The function D{u) has NP part which was measured in and has exponential form, 

/}(«) = 1^(0) exp(-U). (11) 

Finally, the string tension a can be expressed through D{u) ( and higher correlators) 
and at least for static quarks D{u) yields dominant (up to few percent p, ^]) contribution 

to (T, 

a = ^ J D{u)(fu. (12) 

Having in mind all the above relations, one can now average over gluonic fields in (0) to 
obtain 

where the EQL proportional to {{A"-)) is denoted by L^eqli 

LfoL = C / d'xdS (x) ^V(x) V+(y) 'i^as{y){A^£{Mai{y)h'^H:} (14) 



L%L = It / d'xd^d'z (x) ^v(x) ^^+(y) ^i^,M Vi>(^) V/.(^)x 



X (A;jV)^S'(!/)^S'(^)>Ti'j'7i?Tji' (15) 
Average of gluonic fields can be computed using (^,(|3), (||) as 

g'^{A^J^\x)A^^\y)) = ^^^^ [ duia^{u) [ dvka^{v)D{u - v){6^Jik - diuhf,)- (16) 



Jo 

The corresponding expression for the triple average (A^) is given in Appendix 1. 

As it was argued in pGf the dominant contribution at large distances from the static 
antiquark is given by the color-electric fields, therefore we shall write down explicitly 
L^Q^(e/) for this case, i.e. taking /z = z/ = 4 in (0), (|I6|) , while the general case, also for 
generalized gauge (^) is given in Appendix 2. As a result one has 

L%Liel) = ^Jd'^J d'y ^^l(x) / V(x) (y) ^i^ae{yW:h^^^J{x, y) (17) 
where J(x, y) is 

r'x ry 

J{x,y) = / dui / dviD{u — v), i = 1,2,3. (18) 
JO Jo 

Note that we have omitted everywhere for simplicity the upper index (Y) in ^^^^^x), but 
it is implied, since otherwise both (|1^ and (^) are not gauge-invariant. 



3 Bosonization of Effective Quark Lagrangians 

We shall separate out white bilinears in (|1^ following the standard procedure given for 
the general L^eql [H,p7|, 



(n) 

fg 



and introduce the isospin generators t*-"'' 



ni-l 



V " 1 1 

E4",'4f = ^^//„;*'°' = ^i. (20) 



Hence the bilinears in (|T7|) can be written as 



n)-l 

"^Lli^.yWiiy.x) = 2 y: (21) 

n=0 

where we have defined 

V1/W(X,1/)^ fi^tMtf^'i^ae{y). (22) 

Now one can use the Fierz transformation (see Appendix 3 for more details). 

^(4) -,(4) _ i ^231 

Ial3la'l3' — a 2^ ^al3'^a'f3 \^'^) 
^ k=l 

with 

0(^)0(°) = (a,fc)(a,fc) - {aik){cT4k) - (^fc4)(crM) (24) 

Introducing ( PT| ) and ( P3[ ) into (P!7|), one obtains 

4'l?L(e/) = -Jd'xJ d'y^^-''\x,y)¥--'\y,x)J{x,y) (25) 
where we have defined 

^^"-''Kx^y) = l¥:J{x,y)Oi'J, J ^ ^J. (26) 

Bosonization is now done in a standard way using identity (signs and indices of summation 
and integration are suppressed) 

g-M-J* ^ J ^^g^ jf^Dx exp[-xJx + ^'^JX + (27) 

and hence the partition function with the only Gaussian contribution L^^^ assumes the 
form 

Z = J DiPDtfj+DxexpLQML (28) 



where the effective quark-meson Lagrangian is 



'(2) 
^QML 



and the effective quark-mass operator is 



(29) 
(30) 



n,k 



In a similar way one can bosonize the term Leql- However the computations are more 



lengthy and we shall present the result in a separate publication |jT3| 



4 The effective chiral Lagrangian Lecl and stationary 
point analysis 



We are now in position to integrate over quark fields in (^) and obtain Effective Chiral 

(2) 

Lagrangian Lecl- The result is 



with 



L 



(2) 
ECL 



ix) 



Dxe 



d'xd'y Y: X^'^''\x, y)J{x, y)x^"''\y, x)+ 

n,k 



(31) 



+ N^tr \n[{id + im)S{x - y) + iM{x, y)]. (32) 

In ([3^ ) we have taken into account that M is colorless, and the sign tr refers to the 
summation (integration) over Lorentz, flavor indices and space-time coordinates. 

As the next standard step one finds the stationary point equations to determine the 
ground state values for the auxiliary bosonic fields x- Taking functional derivative of 



Lecl with respect to one obtains 



2x^-''\x,y)J{x,y) = -zNMS{x,y))0'^'¥'^^]J{x,y) 



where we have defined as in 14 



S{x, y) = -[{id + im)i+ iM]'].. 



(33) 



(34) 



The set of (nonlinear) equations (^31), ( ^41) is the central result of the present paper. 
In what follows we shall study the properties of the solutions and compare this result 



to the previously obtained equations in []I4[, where another method was used -large Nc 
approximation in the Dyson-Schwinger equations for the heavy-light system with one fixed 
flavour. In this case one should take nj = 1, and = Moreover, only scalar part 
(since it is dominant at large distances from the heavy source [H, 26] was considered, 
hence one can write instead of ([30| ) 
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M ^ Mo(x, y) = -^X{x, y)J{x, y) (35) 
and (|33|) reduces to the equation 

iMo{x,y) = {jiS{x,y)^i)k=iJ{x,y) (36) 
where we have used the relation. 

{l.S^.U = 7S^/3'a'7i'i = \Y. 0^^^tr{SO('^). (37) 

The equation ( ^61) is exactly the same as Eq. ( pTSf ) in [l^ where color-electric field com- 
ponent is retained, namely in the full answer 

iMo{x,y) = J{x,y)-f^S{x,y)'y^- JiklkS{x,y)-ii (36') 

one keeps only /i = 4, = 1. As it is follows from the definition of S", Eq . (|3^ , one has 
another equation to complete a full set 

{-id - im)S{x, y)-i I M{x, z)S{z, y)d^z = 6^'^\x - y) (38) 



Let us look more closely at the scalar part of mass operator, M(j{x,y) Eqs. (|35D, (pq). 
The properties of the kernel J{x,y) ([TB|) have been thoroughly investigated in |T^-|T^ 
and it was shown there, that when x is close to y, then J{x, y) is growing linearly at large 
|x|, e.g. when for simplicity D{u) is taken in the Gaussian form 



D{u) = D{0)e (39) 

then 



J(x ~ y, |x| ^ oo) = |x|2Tgv^D(0)e-("^-^^) (40) 

whereas 

a = ^ y D{u)d\ = 27tT^D{0). (41) 

Moreover S{x, y) (|3^ ) (where only Mq is retained) at large distances displays the properties 
of the smeared 6 -function (see |]l4|-|jl6| for discussion and numerical estimates) 

74^(x,y)74~5('Hx-y) (42) 

and as a result the product J{x,y)'y4^S{x,y)'j4 behaves linearly in |x| at large |x|, in such 
a way that in the equation ( ^8|) one has 

Mo(x, z)S{z, y)d^z ^ a|x|5(x, y). (43) 

Thus for |x| ^ Tg one has 

Mo(x, z) ^ a|x|5(^)(x -z){l + 0{TJ\^\) (44) 
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and 6^^\x — z) is smeared off at distance of tlie order of Tg. To proceed we disregard first 
in the sum (|30D all terms except for the scalar and pseudoscalar fields, 

M{x,y) = xs{x,y)—^ + Xn{x,y)ij5J{x,y) (45) 

with 



The form (45) can be equivalently parametrized in a nonlinear way as follows 

M{x,y) = Ms{x,y)IJ{x,y),U = exp{i^^4>),(i){x,y) = (j)Ax,y)t-l' . (46) 
Now using normalization property 

tr(t(")t(™)) = ^6nm, m = 0, 1, . . . , - 1, (47) 

one easily obtains 

ltr(MM+) = lixl + xl)J\x, y) = M^nf (48) 
and hence the first term in (^) can be written as 

ix1 + xl)Jix,y)=2nfJ-\x,y)Mlix,y) (49) 
and the total Lagrangian (^) is 

L%j^{MsA) = -2nj{J{x,y)Y^Ml{x,y) + N,tr\og[{td + im)i + iMsU]. (50) 
The stationary point equations assume the form 



6L 



(2) 

ECL _ , I T/ 



4nj( J(x, y))-'Ms{x, y) - N,tr{Sie'^^^) = (51) 



SMs{x,y) 

with 

S{x, y) = -[id + im + iMsU]~l, (52) 
5L 



(2) 

-^-^ = NMSMse'^^'^^,). (53) 
b(t)[x,y) 

The solution with = 0, M5 = mP satisfies(|3) while(|l]) may be rewritten in the form 

iMf\x,y) = ^trSJ{x,y) = N,{j,Sj,)J{x,y), (54) 

where {T)sc means the scalar part of operator, as defined in {^7\), and one can see that one 
recovers Eq. (|36|) , derived before in [|14| in a different formalism. 
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5 The derivative expansion of L^ecl 

In this chapter chiral Lagrangian will be written as a series in powers of derivatives of the 
field U{x,y). A similar procedure for the local case and in absence of confinement was 



systematically done in and more recently in [^. For the case of instanton 

model it was done in . In all cases at some moment the local limit of the resulting chiral 
Lagrangian is done, and confining properties of the kernel Ms{x,y) are not taken into 
account. In what follows we keep both (j) and Ms nonlocal and the confining property (^41) 
is exploited. It should be noted that nonlocality of the field is a necessary consequence 
of its quark- ant iquark structure and this structure is lost when localization approximation 
is done. Since the radius of the pion is around 0.6fm, it is only for small momenta that 
localization procedure is justified. We now turn to the second term on the r.h.s. of ( [501 ) 
which contains pionic field and make an expansion of its real part in powers of derivatives 
of the field U. Note that imaginary part of effective chiral Lagrangian was studied in 
0)0) 01 it starts with the terms of the 5-th power in pionic field; it will not be 
studied below. Defining the real part of the pionic effective action i?eLe//[vr], one has 

ReLeffln] = logdet^j^, (55) 

where notations are used 

D = id + im + iMsU,DQ = id + im + iMs. (56) 

Moreover 

D+D = -a^ - d^iMsU^h,, + (m + Msf + mMsiU^ + tj) (57) 
D+D^ = + (m + Msf. (58) 
To simplify we put m = and rewrite (ESl) as 



ReL.jfln] = -^trlog(l - God^MsU^^^), (59) 

where we have defined 

Go(a;,y) = (-92 + M|)-J. (60) 

In (^ODthe sign tr implies the sum over Lorentz and flavour indices and integral over 
coordinates, having in mind that every factor appearing in (|59D under the logarithm sign 
should be considered as a matrix in the space-time coordinates. Expanding logarithm in 
(|59|) and keeping only 2nd and 4th order terms in d^U^ , one obtains 



+ ^God,U^GodMGod^U+GodpU{6,J^p + 6,^6,^ - S^J,^)] + ... (61) 

Here tr implies summing over flavour indices and integration over coordinates, and U = 
Ms{x,y)e^'^^^''^\ We consider now the first term on the r.h.s. of (pl|) and write explicitly 
the coordinate part of tr 

Rel'^^jfl'ir] = N,trfRe[ J d^xGo{x, y)d^U+{y, z)d'^zGo{z, u)d\d^U{u, x)]. (62) 
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In the limit, when Ms{x,y) becomes local, 



U{x, y) ^ ¥{x - y)Ms{x)e^^^''^ = 5\x - y)U{x) 

one has 



ReLlfflTT] = N,trfRe[J d''xd'yGo{x,y)d,U+{y)Go{y,x)d,U{x)]. (63) 



One can see that Eq. (|63D yields a nonlocal chiral Lagrangian, the nonlocality being given 
by the range of the quark Green's function Go{x,y). In case of a model without con- 
finement, e.g.instanton model or NJL model, Go{x,y) would represent the free Green's 
function of a massive quark,with the constituent mass /i created by chiral symmetry 
breaking, /i ^ 0.3GeV. Hence the range of nonlocality in this model case is large, and 
the derivative expansion in powers of dU, which is a standard systematic procedure in 
BjOjIH is justified for small momenta p < fJ,- At this point one should remember that 
in (pTD-(|5^) both Ms{x,y) and U{x,y) (and hence S{x,y), Eq(P^), are defined gauge- 
invariantly with respect to the contour Y, which for simplicity was taken to be X4 axis. 
Therefore the resulting string in M, Ms goes from the points (x, y) (coinciding when Tg 
tends to zero) to the contour Y. Physically it means, that the effective Lagrangian ( |62D 
describes the bosonic field (e.g.pion) in the presence of the heavy quark. It will enable us 
to define in the following chapter matrix elements of heavy-light meson transitions with 
emission of a pion. The case of effective boson Lagrangian for light quarks and antiquarks 
(i.e. without the heavy quark line) will be considered in the second paper of this series. 

6 Pionic transitions in heavy-light mesons 

In this section the main emphasis will be on the pionic part of the quark mass operator, 
which enters the quark-meson Lagrangian as 

AL = t [ d'^xd'^yi/j+{x)M{x,y)ij{y) (64) 



where M{x, y) according to ( ^6]) can be written in the form 

M(x, y) = Ms{x, y)e*^5<^(^'S'). (65) 
In the limit of small Tg one obtains for Ms{x,y) a localized expression 

Ms{x, y) ^ a\-x\5^^\x - y), |x| > Tg (66) 
and the corresponding linearized in Lagrangian is (in the Minkowskian space-time) 

AL(^) = / iP{x)a\x\-f5——i;{x)dtd^x (67) 

where is known to be = 9AMeV. For the pionic transition between heavy-light 
states one should compute matrix elements 

(M2(p2),7r(k)|AL«|Mi(pi)). (68) 
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Neglecting recoil momentum of heavy-light meson one effectively reduces the problem to 
the calculation of the matrix element 



W21 = I i^^{x)'^^,^^^U-^)Sx (69) 



and the decay probability is 

w = 27r\W2i\'5{E, -E2- (70) 

In ( |66D ilj2{^)y4'i{^) quark wave functions of the heavy-light states 2 and 1 respec- 
tively, which can be taken from the solutions of the corresponding Dirac equations, found 
in At this point one can rewrite the matrix element of AL^^^ between two 

stationary quark states iprnyi^n- Using Dirac equations 

[OLp + p{m + Cr\x\) + Vcoul]^n = ^ni^n (71) 

V'm[-Q:p + P{m + a-|x|) + Vcoul] = (^mi'm (72) 

one obtains 

{m\AL^^'^\n) = 77^(m| - 2m75^ + /575(em - e„)^ + 75/3Q;p7r|n). (73) 

In the chiral limit, m = 0, one can rewrite the last two terms inside the brackets in (|73|) 
as 

^^^^d^ ^ ^^~di' ^ ^57M^/^7r. (74) 

One obtains from (|73[) the form of the quark-pion interaction which one usually writes as 
(see and refs. therein) 



AL"^ = c^^tr (^7^75^^ V^), u = -—{ud^,u+ - u+d^u), (75) 

where u = \[tj and g\ is the axial vector coupling of the (constituent) quark. It is easy 
to see that (^) is the first term in the expansion of (^) in powers of the field vr, and it 
is gratifying that the parameter g\ is defined theoretically to be equal exactly to 1 in our 
local approximation,which is in agreement with large argument in pSl . 



7 Concluding remarks 

We have performed a systematic bosonization procedure starting from the QCD La- 
grangian and have derived the nonlocal chiral Lagrangian and its limiting local form 
for quarks interacting with pionic field. By the method of construction, the resulting La- 
grangian is applicable for quark-pion interaction when the quark is coupled by the string 
to the heavy antiquark. Therefore our results can be immediately applied to the pionic 
transitions in the heavy-light mesons. The one-pion transitions have been studied using 
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the form (|7S|) in [Q, and it was found from comparison to experiment that the values of 
g\ around 0.7 are preferred. It makes it reasonable to study the nonlocal version of the 
chiral Lagrangian ([73D, since nonlocality effectively decreases resulting matrix elements, 
possibly explaining the mismatch between our theoretical value g\ = 1 and observed 
value of 0.7. In a similar way one can obtain matrix elements for double pion emission in 
heavy boson transitions, which are of special interest for {B, D), {B, D*) semileptonic de- 
cays. The case of purely light mesons, can be treated in a similar way, but needs another 
string configurations to be taken into account. This is planned to do in the subsequent 
paper of this series. The author is grateful to F. Gross and I.Musatov for discussions and 
to J.Goity for useful criticism and suggestions. This work was supported by DOE con- 
tract DE-AC05-84ER40150 under which SURA operates the Thomas Jefferson National 
Accelerator Facility. 



Appendix 1 

Vacuum averages of Field Correlators 



For the field operator transported to the point Y, 

F^P(x) = ^Y,x)F,,{x)^x,Y) (Al.l) 
so that gauge transformation has the form 

FS:\^)^U^iY)F!,V{x)U{Y), (A1.2) 

the vacuum average of any product can be expressed through Kronecker symbols, e.g. 

= 6adSbcV{x,y) (A1.3) 

where ab, cd are fundamental color indices and V{x, y) can be easily connected to the 
color trace, i.e. to the FC, 

V{x,y) = ^/(tr(F(^)(x)F(^)(y))). (A1.4) 

c 

Similar rules apply to products of any number of they can be deduced from the 
global SU{Nc) rules for the tensors averaged with the Haar measure, namely one can 
make a gauge transformation 

{F^Ax))a>, - n+{Y)UFU^))cMYU (A1.5) 

and average over DQ with the usual Haar measure. In this way one obtains relation (H). 
At this point it is interesting to note that the product of three F's may have another 
representation in the SU{3) group, since there one can use totally antisymmetric tensor 

^abc- 
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E.g. one can write the following equality 

^ace etdf — ^abi^cdSef - ScfSed+ 

+SadiScf5eb - 5cbSef) + 
+ 5af{5cb5ed - 5cd5eb) (Al.6) 

Therefore on can use two forms of writing for the product. In the first case one writes 

'P2{x,y,z)6bcSfcSad (A1.7) 

and finds Vi,V2 by multiplying both sides of (|A1.7|) with the corresponding contribution 
of 5-symbols. In this way one arreves at the Eq. (^) in the main text. If instead one uses 
(|A1.6|) instead of one of combinations, or separates ( |A1.6|) out of Eq. (|), one arrives at 



the white (3g), (3g) combinations. 

Appendix 2 

Vacuum averages of ((^)") 

One can use the generalized contour gauge expression for A^lx). 

dz,-^F,p{z) = J dT^,p{z)F,,{z) (A2.1) 
to represent the average of the product of any number of operators as 

{i^lLi^)aibi---{^lin{Xn))a„b„) — 
dT,,,,,,iz,)... r rfr^„.„p„(^„)((F(f;,(^l))a,.....(FSj^n))a„6^^ (A2.2) 

In a particular case of the modified Fock-Schwinger gauge one has: F = OjXj — > Xj 

dT^^p{z) = a^,{z)dzjp^, (A2.3) 

Appendix 3 

Fierz transformations 



In this appendix the derivation is given of Fierz tranformations for combinations of 7^ 
matrices met in the text above. It is based on the clear presentation done in the book 



29| . Note however that we are always working with the Euclidean 7 matrices, therefore 



some details and coefficients obtained below are different from 129 . 
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We start with the general expansion for any 4x4 matrix 

1 



7 



4 



lA = 1, 

M "Ail' 2j • 



Y^CaIa, A = 1,...16, 
A = 1 

/i = 1,2,3,4; A = 2, 3, 4, 5 

A = 6,7,8,9,10,11 
A = 12,13,14,15 
A =16 



(A3.1) 



(A3.2) 



7a = 757m, 
7a = 75 

In (A3.1) one can derive the general representation for any matrices Fmk,Gii, indeed 
from relation 



changing m ^ m', I ^ V and multiplying with Fmm'Gin one obtains 

1 



PmkGil 



4 



J2 Aa{F^aGU{ia: 



ik • 



(A3.3) 



For Euclidean matrices 'Ja one easily obtains = 
= 1 otherwise. Taking F = G = '-f^, one obtains 

iLklti = 7llAA(7V7^)mz(7A)ifc 



1, for A = 12, 13, 14, 15 and 



Ak\/4:k\ 



-{a''){a'') + (7V)(r7') - - {llill} 

A ^rnl ^ik ■ 
^ k=l 



(A3.4) 



The last relation coincides with p3| ) and operators O^'^^ are given in (p4D. Note that 
in the curly brackets in ( |A3.4| ) each product of 7 matrices, (77) (77), has the same order 
of indices as in ol^}ol'l\ In a similar way one can represent the combination of spacial 7 
matrices, n = 1, 2, 3, no summation over n, 

llkla = \{{m) + (7')(7') - (7™)(7'") + (7")(7") + 



+ [an{crnn^,. - {cr^ni^^n ' - (7'7'^)(7'7'^)M^n+ 

+ (7V)(7V)-(7')(7')W^. 
Summing (|A3.5|) over n and adding ( |A3.4|) one obtains 

(7^)^^.(7^). = {(1)(1) - ^(7'^)(7^) - ^(7V)(7V) - (7')(7') W- 



(A3.5) 



(A3.6) 



Now in case of generalized contour gauge as in Appendix 2 one has to make Fierz trans- 
formation of the combination 



{l^mkiinu = jT.^A{l'lAY)ml{lA)^k. 
^ A 



(A3.7) 



Note that scalar and pseudoscalar combinations occur on the l.h.s. of ( |A3.71 ) only for 
/i = z/. 
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